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1. Introduction 

A central issue in modern physics is the investigation of noise as it drastically affects the 
evolution of quantum systems. Such a general phenomenon is of importance in quantum 
information science and beyond, as it addresses a fundamental issue in quantum theory 
[U E]. So far, two main dynamical regimes, Markovian and non-Markovian, can usually 
be distinguished according to the noise time scale (respectively shorter or longer than 
that of system dynamics). Here we address the non-Markovian case. In facts, despite 
Markovian evolutions have been exstensively investigated, in general real noisy dynamics 
are far from being Markovian. Although gaining much interest in the last years both 
in theory and experiment [31 H], a general theory for non-Markovian dynamics is still 
missing. In this dynamical regime, exact master equations have been derived for 
the evolution of a Brownian particle linearly coupled to a harmonic oscillator bath, 
for instance via path integral methods O |6] or phase-space and Wigner function 
computations [3 [S]. Analogous results are derived by means of quantum trajectories, 
either exactly or in weak coupling approximation [9l [101 HB 112] • In the framework of path 
integral methods, master equations have been derived both for initially correlated states 
[131 [H], and for factorized initial states in the case of weak non linear interactions [15] . 
In general, all these master equations cover only few cases and are not simple to solve. 
Indeed, it would be highly desirable to find an approximation scheme fully capturing 
non-Markovian features. In general, different approximations (for example on system- 
bath interaction strength) may lead to irreconcilable dynamics [T]. 
For Markovian dissipation it has been shown that, by exploiting symplectic tomography 
[IB] and using Gaussian probes, unknown parameters governing the dynamics of the 
system can be reconstructed trough a limited number of tomographic measurements [IT] . 
In this paper we extend this tomography-based approach to time-dependent dissipative 
dynamics. We present an experimentally feasible procedure that exploits symplectic 
tomography and which, under the assumption of Gaussian noise, allows to reconstruct 
unknown time- independent parameters that characterize the time-dependent coefficients 
of the master equation. Even though the assumption of Gaussian noise might be seen as 
an idealization, it is actually well fitted for a significant number of models [H |2]. Also, 
small deviations from Gaussianity would introduce small and controllable errors. We 
also note that Gaussian probes are quite straightforward to produce [T8] . 
The paper is organized as follows. In section |2] we briefiy recall the procedure that 
allows to reconstruct cumulants of a Gaussian state (our probe) undergoing a dissipative 
dynamics. In section |3] we introduce the class of non-Markovian master equations we 
will study, and derive the expressions for the first and second time-dependent momenta 
(cumulants). In section |4] we devise two alternative approaches based on quantum 
tomography that allow to reconstruct unknown time-independent master equation 
parameters. In section [5] we apply these procedures to a benchmark model and compute 
the amount of measurements needed. In section [6] we summarize and discuss our results. 
Appendix A contains more details on the reconstruction of the cumulants of a Gaussian 
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state through symplectic tomography. 

2. From tomograms to cumulants 

In [TTj we have introduced a procedure that allows to reconstruct, via a limited 
number of measurements, the time-independent master equation coefficients governing 
the dynamical dissipative evolution of a quantum system. In the Markovian case, it 
has been shown that, using a Gaussian probe, the required number of measurements 
is at most ten. In facts, the dynamical evolution of a Gaussian state is completely 
determined by the evolution of its ffist and second order cumulants, which are 
measurable quantities. The unknown master equation coefficients enter the dynamical 
equations of the cumulants, hence can be retrieved by simple inversion, once the latter 
are measured. The cumulants can be obtained by using symplectic tomography. Indeed, 
given the Wigner function of a Gaussian state, it can be measured along lines in phase 
space (i.e. by performing its tomographic map). This allows to relate the cumulants 
to points on the tomogram. By choosing the lines in phase space corresponding, 
respectively, to position and momentum probability distributions, one needs at most four 
points along each tomogram (i.e. line) to retrieve the ffist and second cumulant of the 
associated variable. The same procedure, applied on a line inclined by tt/4 in phase space 
allows to retrieve the covariance of the two variables by measuring at most two points. 
Hence, given a Gaussian Shape Preserving (GSP) Markovian master equation, by using 
Gaussian probes one can retrieve at any time its evolved cumulants via a finite amount 
of measurements. In the following, we will refer to this as the tomograms-cumulants 



(T-C) procedure. More details of the T-C procedure are provided in Appendix A As the 
key ingredient of the T-C procedure is the preservation of Gaussianity, it is therefore 
a natural step to investigate how this procedure can be generalized and extended to 
more involved non-Markovian scenarios, still preserving Gaussianity. Indeed our aim 
is further supported by recent work fi9[ [20l [2T] . in which it has been proved that it 
is in principle possible to make tomographic measurements of the probability densities 
associated to every quadrature in phase space (for example in quantum optics it could 
be realized by means of homodyne detection). As a final remark we note that other 
methods to measure the covariance matrix of Gaussian states have been discussed in 
[22] . However in this case the amount of required measurements is much higher, the 
focus being on the reduction of experimental errors. 

3. Non-Markovian master equation 

We will focus on the class of master equations of the form [T] : 



dt 



mm)] = -^ [H,m] + ^E ([^^w/^W'^^w] + [m,pmHt) 
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where the generator i^(t)[-] depends on time. In particular we shall investigate the 
following class of time-dependent master equations 

^ = 4 [i/„. P«)] - '-^^^^ iimp^pm]^ '-^^^ \P, pm+ m)] 

- %^ [9. I?.p11 - %^ [p. [p.^tJllH- ^ (I?. Ip.pWD + |p. B.^WU) -(2) 

The master equation (|2]) is obtained from the general form (fTl) by choosing a system 
Hamiltonian of the form 

H = Ho + - (qp + pq) , i^o = ^p' + ^-g', (3) 

and linear Lindblad operators Vj{t) 

V,{t) = a,{t)p + b,{t)q, J = 1,2. (4) 

Using equations ^ and Q in ([I]), the time-dependent master equation coefficients in 
Q hence read 

i=l,2 j=l,2 

j=l,2 j=l,2 

where the Dij{t)'s and A(t) are, respectively, the so-called diffusion and friction 
coefficients. The master equation in ^ represents a natural generalization of the time- 
independent master equation introduced in [23] which describes a GSP evolution of 
a quantum state. In the time-dependent case a wide range of models obeys a GSP 
master equation of the form (|2| [6l [3 [U [101 HB 1121 HSl [Hj. As a side remark we 
note here that the literature about non-Markovian master equations may lead to some 
ambiguity. Indeed some authors classify as non-Markovian only those master equations 
whose generator contains a convolution integral. It has recently been proved in [21] 
that these generators can be mapped into convolutionless ones, following a so-called 
local approach. Non-Markovianity becomes then characterized by the dependence of 
the convolutionless generator on t — to where to is the initial time. According to this 
approach a time- dependent convolutionless generator as equation ^ could be considered 
Markovian. However, following a consistent part of literature, e.g. [H [3, [HI dUl IHl 1151 ES] , 
we term non-Markovian convolutionless time-dependent generators as the one in (|2|. 

3.1. Evolution of the cumulants 

We begin by writing the dynamical equations, obtained from (|2]), for the cumulants of 
a Gaussian state. They can be expressed in compact matrix form as 

j^S{t) = {M-\{t)h)S{t\ (6) 

jX{t) = {R-2\{t)h)X{t) + D{t), (7) 
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where /2(3) is the 2(3)-diinensional identity matrix. The vectors S{t) and X{t) 
correspond, respectively, to the first and second order cumulants 



S{t) 



^/muj{q)t 



Vn 



{p)t 



X{t) 



n 



\ 



/S.p1/{muj) 



(8) 



the matrices M and R contain the Hamiltonian parameters 

f 26 2a; \ 



M 



5 CO 
-u —6 



R 



V 



-25 



-2uj 




(9) 



and, finally, D{t) is the diffusion vector 

/ muDgg ^ 

m = -^ t ■ (10) 

As in the time-independent case, the first cumulant dynamical evolution ^ depends 
only on A(t) also entering, together with the diffusion coefficients, the second order 
cumulants equation ([T]). Hence, the friction coefficient can be retrieved by inverting ([6]). 
By carrying out a formal integration one obtains the following expression 

'^.(0)' 



dt'X(t') = In 

\S,{t) 

where the suffix j = 1,2 labels the two components of the vector 

S{t) = e-'^^'Sit). 

Analogously, ([T]) can be rewritten as 

|a-W = Dit). 



where 



X{t) = e2/o'^*'Mt')e-*«x(t), D{t) = e^/o '^*'M*')e-*«Z}(t). 



(11) 

(12) 

(13) 
(14) 



It can be shown [23] that the transformations in (12) and (14) are always invertible. 



provided one sets the quantity r] = y/S"^ — u"^ to iQ whenever 77^ < 0. The formal 

(15) 



solution of ( 13 ) is given by 



X{t) = X(0) + / dt'D{t'). 



Inverting the transformation in (14), one can write 



rdt'e-2//"^*"^(*")e(*-*')^D(t') = X(t) - e*^e-2/o«'*"A(*")x(0). 
Jo 



(16) 



Both in (11) and (16), measurable quantities appear on the right hand side, whereas 



the (unknown) coefficients are on the left. 
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4. Reconstruction of time independent parameters 

Hereafter we will assume that the master equation coefficients (MECs) A(t), Dgq{t), Dpp{t) 
and Dqp{t) have a known functional form. This implies that the non-Markovian mas- 
ter equation with certain expressions for the MECs has been previously derived within 
some approximation scheme (e.g. by means of a microscopic derivation and subse- 
quent dynamical assumptions). The time-dependent MECs are thus function of a set of 
time-independent parameters (TIPs) whose value is a priori unknown. In this section 
we propose two alternative procedures aiming at reconstructing the TIPs by means of 
symplectic tomography. 

4.I. Integral approach 

Here we introduce an approach based on formal integration of the dynamical equations 



(11) and (16). The right hand side of both equations involves experimental inputs 
{S{t) and X{t)) and known Hamiltonian parameters {R). The left hand sides, once the 
MECs are known, can be regarded as functions of the TIPs. Hence the set of TIPs 
can be in principle obtained by inverting these relations. Unfortunately, in general, an 



analytical inversion of (11) and (16) may represent a highly involved task. In facts, even 



if we assume a known MECs time-dependence, we could be unable to either compute 



analytically the integrals on the right hand side of (11) and (16), or to invert the 
equations or even both. All these problems can be anyway overcome by resorting to 
numerical computation. To provide an example of how to apply this procedure, in the 
next section we will apply it to a specific model. 

4-2. Differential approach 

The tomographic T-C procedure recalled in section |2] allows us to measure not only the 
cumulants of a given Gaussian state but also their first time derivatives. Indeed, we can 
estimate the derivative through the incremental ratio by measuring each cumulant at 
two different times t and t + 6t. For instance 

i^A,? ~ MiizM, (17) 

dt ^ 5t ^ ^ 

where the amount of time 6t is defined as the smallest time interval which allows 

to experimentally distinguish two different values of the given cumulant. Once we 

substitute derivatives with their approximations, the two sets of equations ^ and (JTJ) 

are not differential anymore. Being the cumulants and their approximate derivatives at 

given times experimental inputs, we insert in ^ and ([T]) the time-dependent coefficients 

A(t) and D{t) (and in case m{t), u}{t), S{t)), which involve the unknown TIPs and in 

this way the two sets of equations ^ and ([T]) reduce to algebraic equations which can 

always be solved, at least numerically. Due to this simplification, this approach is also 

suitable when dealing with more complicated generators than the one considered here 

in equation ([2]), e.g. when the Lamb shift contribution is explicitly taken into account 
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such that the Hamiltonian parameters become time dependent: m{t), u}{t), 5{t). We 
note that the two sets ^ and ([T]) consist of five equations which must be fulfilled at 
any chosen time. Therefore, considering them at different times, we can derive a system 
made up of an arbitrary number of equations. The number of equations must then be 
chosen as the minimum amount of equations needed to uniquely determine the TIPs, 
which is clearly model dependent. In the following section, we show how to apply this 
procedure to a benchmark model. 

5. A benchmark model 

In the following we will refer to a specific model of a quantum Brownian particle 
discussed in [25]. We will show how to apply in this specific case the two general 
procedures presented in the previous section. The model consists of an Ohmic reservoir 
made of harmonic oscillators, linearly coupled to a single harmonic oscillator of frequency 
u (our system particle) through the coupling constant a, with a Lorentz-Drude cut-off 
[1] Uc and at temperature T. Starting from a superoperatorial version of the Hu-Paz- 
Zhang master equation [5] , a secular master equation of the form ^ is obtained in the 
weak-coupling limit (up to the second order in a), with the following coefficients: 

A -n n -n rnujDgq _ Dyy _ A(t) 

- u, i^,p - u, ^ _ ^^^ _ ^ , 



2 2 

A(t) = 4^^ { 1 - e-"^* cos(wt) + — sin(u;t) 






UJc 



cos[ujt) sin[ujt) 



the last having this form for high temperatures T. This master equation is of Lindblad- 
type when the coefficients A(t) ± A(t) are positive at all times. The Lindblad-non 
Lindblad border as a function of the temperature T and the frequency cutoff Uc has 
been analyzed in |25]. As Gaussianity is preserved, by choosing a Brownian particle 
initially in a Gaussian state, the T-C procedure can be employed at any time. The 
coefficients X(t) and A(t) reach stationary values for t ^ l/coc 

a'^u'^u ^ , , 2a'^u'^ kT .. 

ojj + 00^ u^ + u^ n 

In this specific model, the unknown TIPs are the coupling constant a, the temperature 
T and the frequency cut-off Uc- Usually, when studying quantum Brownian motion, one 
assumes Uc/u) ^ 1, corresponding to a Markovian reservoir, with Uc — ?■ C)0. In this limit, 
the thermalization time [25] is inversely proportional to the coupling strength, while for 
an out-of- resonance engineered reservoir with Uc/u) ^ 1 (i.e. highly non Markovian), 
the thermalization process is slowed down. 

5.1. Example: integral approach 



Here we apply the integral procedure (section 4.1 ) to the benchmark model. In this case 



the left hand side of (16) is not simply analytically computable. Thus we must use (11) 
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a 



a 



0.012 



0.01 



0.008 




cojoj 



(a) 




cojoj 



Figure 1. We show how to indirectly measure the time independent parameters tUc 



and a , in two different regimes respectively close to the Markovian (figure 1(a)) and 



non-Markovian dynamics (figure 1(b)). Each line refers to an experimental measure 
of In (Sj{0)/Sj{t)] at a specific time ojt. In both regimes, the time independent 
parameter values are found at the intersection of the two lines. In the example close 
to the Markovian regime (figure 1(a) I if we measure 3.03 • 10"'^ at tuti = 0.5 (solid 
line) and 9.70 • 10~^ at ujt2 = 10 (dashed line), we retrieve a^ — 0.01 and Wc/w ~ 10. 



Analougsly, for the example close to the non Markovian regime (figure 1(b) I if we 



measure 4.55 -10 at ujti = 0.5 (solid line) and 1.84 -10 at a;i2 = 10 (dashed line), 
we retrieve a^ = 0.01 and cjc/w = 0.5. 



to reconstruct all the TlPs it involves, and then numerically integrate the left hand side 



of (16). 



The left hand side of ( 11 ) is given by 



dt'X{t') 



Q?lJllJ^ 



{ul + a;2) 



2^2 



Lot 



UOI + ^2 



UJ' 



.Wn 



UJ 



-iOct 



1 — cos(a;t) H ^ — sin(u;t) 



(20) 



By using (11) and (20) we obtain the following trascendental equation for the coupling 
strength a 



,2\2 



a 



In 



S,{t) 



UJ^Up' 



ut^^±:^ - 2^ 



+e 



-UJct 



2^cos(a;t) + ^ 

UJ 



UJ^ 



UJ 



UJ 



UJ^ 



sin (cut) 



(21) 



where the ratio Sj{0)/Sj{t) is the experimentally measurable quantity. Hence, by 



performing two distinct measurements of this ratio we can evaluate (21 ) at two different 



times. We thus obtain a system of two numerically solvable equations, which allows us 
to retrieve the time-independent parameters a and Uc- To provide a concrete evidence 
of the validity of this procedure, we show two numerical examples in figure [1} Indeed 
we retrieve the TIP a^ = 0.01 in two different dynamical regimes, respectively close to 



the Markovian (figure 1(a)) and non-Markovian (figure 1(b)) limit. The last missing 
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9 



parameter is the temperature T entering the coefficient A(t). By using (16), we obtain 
the following equation: 



kT 






X 1 



-LOct 



cos{ijjt') 



u 



UOr. 



sin(a;t') 



-1 



(22) 



where j = 1,2,3 denotes the vector components, and (e*^* ^^^)i ^^^ ^^^ matrix 
elements of the matrix e*-*"* ^^. The exphcit expression of the integral appearing on 



the first hne of (22) is provided in (20). In general, the remaining integrals are 



not analytically computable. However, since all the parameters involved have been 
previously reconstructed, these integrals can be computed numerically. 



Number of tomographic measurements Let us now explicitly compute the number 
of tomograms needed to apply the integral approach to the benchmark model. To 
reconstruct a and Uc each of the quantities Si^2{t) = {^~^'^^S)i2 (^) i^ust be measured 



once but not at the same time, as shown by (21) and figure LLJ Each Si^2{t) is a 



function of the first cumulants of both position and momentum. According to the 



T-C procedure (see Appendix A), the reconstruction of a first cumulant involves at 
most four tomographic points. Thus a and Uc can be obtained via, in the worst case, 
sixteen measures. Furthermore, being e~*^^ an orthogonal transformation, e~^^^ Sj{t) is 
by itself a first cumulant along a time-dependent direction in phase-space. Hence, if 
time-dependent tomographic measurements (i.e. measurements in a frame rotating as 
e~*^^) are allowed, the number of required tomograms decreases to eight, as we would 
only need a single first cumulant {Si(t) or 5*2 (t)). 

To measure T we should evaluate one of the second cumulants at a given time. Following 
the T-C procedure this amounts to two tomographic points. However, the required 
second cumulant has been already obtained when reconstructing the corresponding first 
cumulant, hence the temperature can be retrieved without further effort. This argument 



also holds for time-dependent measurements. In facts, the reconstruction formula (22) 



has been derived from ( 16 ), which can be recast in terms of the variances in the rotating 



frame. The temperature can be then obtained using the variance along the same time- 
dependent direction of the measured ffist cumulant. 

In conclusion, in order to implement the integral approach in the benchmark model, 
according to whether we can perform time-dependent measurements or not, we need 
eight or at most sixteen tomographic points. 
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a 



cojoj 




cojoj 



Figure 2. As for the integral approach! we obtain an indirect measure of tlic time- 
independent parameters Wc and a^, both in the almost Markovian (figure 2(a) ) and in 



the almost non-Markovian regime (figure |2(b)| . Each line refers to an experimental 



measure of 



^{Q)i 



( 



{P)t 



{q)t-{q)t+St 



St 



at a specific time ujt. In the almost Markovian 
example shown in figure 2(a) if we obtain 9.52 • 10"'^ at tuti = 0.5 (solid line) and 
9.90-10"'^ at ujt2 — 10 (dashed line), we retrieve a^ = 0.01 andwc/w = 10. Analogously, 
for the almost non-Markovian case in figure 2(b) if we measure 2.59 • 10"'* at ujti = 0.5 
(solid line) and 2.01 • lO"'^ at Lut2 = 10 (dashed line), we retrieve a^ — 0.01 and 

LOc/iO — 0.5. 



5.2. Example: differential approach 



We now skip to the differential procedure (section 4.2). In (18) the dependence on a^ 
and cOc is factorized, hence using ^ one gets 



a^ ~ 



1 f{p)t {q)t- 



It+St 



m 



St 



UjJ.+U^l 



UJt. 



u 



.—LOct 



cos((X't) -I — - sin((X't) 



.(23) 



{q)t 

Since a^ and Uc are time-independent, they can be determined by solving (23) for two 



different times ti and ^2, and looking at the intersection of the two different solutions. 
This procedure requires to measure the cumulants (g)t, {q)t+5t and {p)t at t = ti,t2 and 



to solve (23) numerically, as shown in figure 2 In other words, the first two TIPs 



a 



and cUc, can be determined by measuring six quantities. As for the integral procedure, 
we retrieve the TIP a^ = 0.01 in two cases, corresponding to the extreme dynamical 



regimes, Markovian cUc/uj ^ 1 in figure 2(a), and the highly non-Markovian, Uc/ui <^ 1, 



in figure 2(b) Again, we are left with determining the temperature T appearing in (18). 



To this end, we consider one of the equations of system (16), e.g. that for {2/h)mujDqq{t) 



which according to (17) we reformulate as 

huj h \ 5t 



^ + 2Xmq^ 



^(^{<l,p)t 



2a^ul 



1-e 



-UJct 



u 



cos{ut) sin(ci;t) 



Ur 



(24) 



Equation (24) allows to retrieve T once the cumulants Ag^, ^q^^s^ and (j{q,p)t are 



measured and both a"^ and Uc are known from the previous steps. As an example, at time 
ut = 1 the same value of temperature, kT = lOfiw may correspond to different measured 
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values of the quantity 1/h [m ((Ag^_^^^ — Aq^)/St + 2\{t)Aqf) — 2a{q,p)t] , according to 
different dynamical regimes. For instance, by setting Uc/uj = 10 close to the Markovian 
regime one would measure 0.198 as output whereas, setting Uc/oo = 0.5 close to the 
non-Markovian case, the output corresponding to the same T is 0.067. 

Number of tomographic measurements We now count the number of tomographic 
measurements to apply the differential approach to this example. The reconstruction of 
both a and Uc is based on (23), which must then be evaluated at two different times ti 



and t2, see figure 2 Each evaluation of (23) requires two measurements of the average 



position, at times tj and tj + 6t, and one of the average momentum at time tj, where 
i = 1,2. This implies reconstructing six first cumulants. As each first cumulant requires 



four tomographic points (see Appendix A), the total amount of needed tomographic 
points amounts to twenty-four. 



The reconstruction of temperature T is based on equation (24). We need the variance 



of the position at time t and at time t + 6t, and the covariance of q and p at time t. 



However, according to the T-C procedure (Appendix A), the variance is required to 
obtain the position average. This implies that the second cumulant has been already 
measured during the previous reconstruction, and there is no need to measure it again. 
Retrieving the covariance requires two more tomographic points. 

In conclusion, in order to implement the differential approach in the benchmark model, 
we need twenty-six tomograms. 

Comparison between the two approaches Let us now briefly compare the two 
procedures described in this section. On one hand the differential approach requires 
more experimental measurements compared to the integral one while, on the other, the 
latter procedure is more involved from a computational point of view. Indeed, it may 
happen that to compute the first members of Eqs. (11)-(16) some numerical or analytical 
approximations are needed, thus reducing the accuracy of the reconstruction. In this 
case, the differential approach should be preferred as it is very simple from the point 
of view of analytical computation. Clearly, if the computation of the integral functions 
in Eqs. (11)-(16) does not present remarkable difficulties, the integral procedure proves 
better as, requiring less measurements, it involves a lower number of interactions with 
the physical system. For example, in our benchmark model, the differential approach 
requires twenty-six measurements, while the integral approach requires sixteen time- 
independent measurements or only eight time-dependent measurements. 
One could summarize by saying that the integral procedure is more advantageous in 
terms of number of measurements, but requires the ability of solving potentially involved 
analytical expressions. The differential approach, instead, is more advantageous from the 
point of view of versatility, as it allows to deal in a straightforward way with complex 
generators, at the expenses of a higher number of measurements. In conclusion, the 
choice between the two strategies introduced in this paper strictly depends on the specific 
model under investigation. 
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6. Conclusions 

In this paper we liave proposed an experimentally feasible procedure to reconstruct 
unknown time-independent master equation parameters in a non-Markovian scenario. 
To this end we have adopted an approach based on symplectic tomography. While 
previously the case of Markovian dynamics has been investigated [IT] , here the procedure 
is generalized and extended to the more involved convolutionless non-Markovian case. 
In particular, our analysis is focused on the class of Gaussian Shape Preserving 
master equations. We have addressed the situation in which the time-dependent 
master equation coefficients are analytic functions of some unknown time-independent 
parameters. The key point of our approach lies in using Gaussian states as probes, 
as information on the dissipative dynamics can be inferred via a limited number of 
tomograms. We have proposed two alternative procedures, integral and differential, to 
reconstruct the unknown quantities. In order to provide an explicit example of how 
these different approaches work, we have applied them to a benchmark model made 
up of a harmonic oscillator coupled to a bosonic bath, whose unknown parameters are 
the coupling constant, the temperature and the bath frequency cut-off. In this case the 
number of needed tomograms ranges between sixteen (at most) in the integral approach 
and twenty-six in the differential approach. 

Besides measuring unknown parameters, our procedure proves useful also in case those 
are already known. Indeed it could be employed as a preliminary consistency test for 
the adopted master equation. In facts, the reconstruction procedure assumes that the 
time- dependent master equation coefficients are previously known functions of a set of 
time-independent quantities. This is for example the case of a microscopical derivation 
(and related approximations) of the master equation. In this perspective, the agreement 
between the measured and theoretically expected time-independent parameters provides 
a necessary validity condition for the adopted approximation scheme. Along the 
same line of thought one could also envisage an extension of this approach to the 
reconstruction of the whole set of time-dependent master equation coefficients. This 
could provide a sound, reliable and complete experimental check of the goodness of the 
approximation scheme underlying a master equation. A detailed investigation of this 
wider scenario will be subject of a distinct study [26]. Our proposal opens up several 
interesting questions which are going to be the subject of further future investigation. In 
facts how our approach can be recast within an estimation theory perspective represents 
a relevant open scenario. A similar analysis has indeed been performed in |27j for 
a single-parameter Markovian master equation. Another relevant point to investigate 
is whether the proposed protocol can be enhanced by employing entangled Gaussian 
states as a probe. Finally, whether or not the proposed procedure can be generalized and 
employed in presence of memory kernels is a challenging question. Indeed, reconstructing 
the unknown parameters of Gaussian noisy evolutions with memory represents both a 
highly involved and interesting task. 
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Appendix A. The T-C procedure 

Tomographic maps [El [28] allow to reconstruct the state or some other properties of a 
physical system, both in a classical and in a quantum regime. In general, tomography- 
based techniques stem from a probabilistic perspective. Indeed, given a quantum state 
p{t) its Wigner function [291 EO] provides a generalization on phase space of a classical 
probability distribution and is defined as 

W{q,p,t) = — dyexpl^^jp{q-y,q + y,t). (A.l) 

The previous equation can be read as a map between real phase-space functions and 
density matrices. In particular, whenever the dynamics of a quantum system initially 
in a Gaussian state obeys the master equation ^ the associated Wigner function is a 
Gaussian function itself and reads 



1 



W{q,P,t)= , exp 

2ny/AqiApf-a{q,p)f 



It 



^Qtip - {p)tr + ^Ptii - {q) 



2[Ag2Ap?-a(g,p)?] 



2a{q,p)t{q- {q)t){p - {p)t) 



(A.2) 



2[Ag2Ap2-a(g,p)2] 

Given the Wigner distribution of a quantum system, the Radon transform [31] 
represents the key ingredient to perform a tomographic analysis. This invertible integral 
transformation allows to retrieve the marginal probability densities of the system, i.e. 
the probability density along straight lines in phase space: 

X - fiq - up = 0. (A.3) 

The formal expression of the Radon transform, for a generic quantum state, is then 
given by 

w{X, /i, v) = {6 {X — pq — up)) = / W{q,p, t)5 (A — /ig — up) dqdp. 

(A.4) 



In case of the Gaussian function (|A.2|) it becomes: 

\2 



wiX, u, u) = ——=. — exp 



(A - /.(g), - u{p)tr 



2[AgiV + Ap?z/2 + 2a(g,p)tH 

(A.5) 

The second cumulants always obey the constrain Ag^/i^ + Ap^z/^ + 2cr(g,p)t/iz/ > as 
a consequence of the Robert son- Schrodinger relation |32j. This matrix inequality is a 
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generalization of the Heisenberg principle. The advantage of using Gaussian probes 
to investigate a dissipative dynamics arises also within a statistical perspective. One 
could in fact wonder whether, due to experimental errors, a violation of the uncertainty 
principle might be observed. This may happen if measurements are performed on 
states almost saturating the Robertson-Schrodinger inequality, i.e. on the minimum 
uncertainty states which are pure. However, our measurements are performed on states 
undergoing a dissipative evolution i.e. on states typically far from being pure hence 
from saturating the inequality. Furthermore any additional noise of statistical origin 
will have the effect of moving the reconstructed state further away from the boundary, 
as noted in [22]. 

Let us now consider the tomograms corresponding to the position and momentum 
probability distribution functions {X = q and X = p): 



zu{X, 1,0) 



n7(X,0, 1) 



1 



AqtV2^ 

1 
AptV2^ 



exp 



exp 



(x- 


- m' 


2Agf 
(X - {p)tf 



2Ap? 



(A.6) 



(A.7) 



The lines individuated by the choices (/i, z/) = (1,0) and (/x, z^) = (0,1) correspond to 
tomograms depending on the time average and variance respectively of position and 



momentum. In order to determine the latter quantities we must invert (A.6) and (A.7) 



for different values of X, i.e. for a given number of points to measure along a tomogram. 
By considering first the direction fi = 1, u = 0, and by inverting (A.6), we obtain: 

1 



(X-(g),)^ = 2Ag,^ln 



(A.8) 



w(X, 1, 0) AgiV27r 

t then we need only the value of the tomogram zu{0, 1, 0) to 

t as a function of Ag^: 



If we know the sign of 

get {q)t, otherwise we need another point. In this way we get 



±Aqtj2\n 



I 



w{0,l,0)AqtV2^ 



(A.9) 



Using (A.9), (A.6) becomes an equation for Ag^ only, and it can be rewritten as 

1 



2Aqt In 



Ct7(X, 1, 0) AgtV2^ 



XT^qtx 2\n 



w{0,l,0)AqtV2^ 



(A.IO) 



This equation is trascendental, therefore we can only solve it numerically. For each X 
and corresponding w{X, 1, 0) there may be two values of Aqt satisfying the previous 
equation. In order to identify one of the two solutions, it is enough to consider two 
points, {{Xi, 'G7(Xi, 1, 0))} and {{X2, tu{X2, 1, 0))}, and to choose the common solution 
for the variance. Hence, whether we know or not the sign of the average {q)t, we need 



three or four points to determine {q)t and Aqt in (A.6). Analogously, we need other 
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three or four points for {p)t and Apt in (A. 7). We now compute the covariance a{q,p)t- 
To this end we consider the tomogram: 

" / \ 2 

zu [ X, —= , —pz = exp 



^ V2 ) 



Aq^ + Ap^t+2a{q,p)t 



.(A.ll) 



This is a Gaussian whose average is aheady determined. Indeed, according to the 
previous steps, we need two more points of this tomograms to determine the spread 
{Aq^ + Apf)/2 + cr{q,p)t, from which we can retrieve (7{q,p)t- 

Hence, the first and second order momenta of a Gaussian state can be measured 
at an arbitrary time t by means of eight or at most ten points belonging to three 
tomograms. 
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